Breakdown of the Luttinger sum rule within the Mott- Hubbard insulator 
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The validity of the Luttinger sum rule is investigated within the prototype tight-binding model of interacting 
fermions in one dimension, i.e., the t-V model including the next-nearest neighbor hopping t' in order to break 
the particle-hole symmetry. Scaling analysis of finite-system results at half-filling reveals evident breakdown of 
the sum rule in the regime of large gap at V 3> t, while the sum rule appears to recover together with vanishing 
of the Mott-Hubbard gap. 
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The Luttinger theorem |ll|,|2t] is the essential building block 
supporting the concept of the Fermi liquid (FL) as formulated 
by Landau yj]. The Luttinger sum rule (LSR) in a homo- 
geneous system relates the Fermi volume to the density of 
fermions irrespective of the presence of the electron-electron 
interactions. Stimulated by experiments on novel electronic 
materials with strongly correlated electrons which indicate 
possible deviations from FL scenario and from LSR fl, the- 
oretical studies of the validity of the LSR and its limitations 
have intensified. 

In a metal at T = the Fermi surface is located by the poles 
of the Green's function (GF) G(k, oj = 0) and the LSR shows 
that the Fermi volume being equal to the density of electrons n 
is unchanged by the interaction. It has been pointed out [5, 6] 
that the original derivation [ 1 ] and the LSR can be general- 
ized as well to insulators where the corresponding 'Luttinger 
surface' (LS) is defined by zeros G(k, 0) = 0. Such LSR con- 
cept becomes of interest, but as well easier to test in strongly 
interacting electrons and in Mott-Hubbard (MH) insulators in 
particular. It has been recently applied to spin ladders On 
the other hand, there are several indications that LSR might be 
violated within the MH insulators in general JH 0, OH . The 
argument is based on the observation, that the LSR is satisfied 
only for a particular value of chemical potential /i within the 
MH gap. It has been shown that for models with the particle- 
hole (p-h) symmetry the latter is the case and the LSR is ful- 
filled [81]. At the same time, it has been realized that the LSR 
should as well apply to finite systems 111 ll Il2ll . This allows 



us to test validity of LSR in nontrivial models of correlated 
electrons lfl3ll . Based on analytical expansion for U/t 1 it 
has been shown on small systems that within Hubbard model 
on a planar triangular lattice (without the p-h symmetry) LSR 
is indeed violated for a range of parameters 111 311 . 

In this Letter we present results of the numerical study 
within the prototype model of interacting fermions in ID, i.e., 
the generalized t-V model. The advantage of such a ID model 
is that it allows for the finite-size scaling to the thermody- 
namic limit. Our results show clear violation of the LSR for 
V ^> t within MH insulating phase which appears to persists 
down to critical V > V c where the MH gap opens. 



In the following we study the extended t-V model 

H = - i I>!+iC J +h.c.)-^(cl +2 c 4 +h.c.) + 

i i 

+ Vy]nini+i, (1) 

i 

where t and t' are nearest-neighbor (n.n.) and next-nearest- 
neighbor (n.n.n.) hopping, respectively, and V is the n.n. re- 
pulsive interaction between fermions. We are interested in 
the MH insulator state which appears at half-filling with the 
electron density n = 1/2. It is well known that the model 
with t' = is equivalent to the anisotropic Heisenberg model, 
which can be solved exactly via Bethe ansatz 1 14]. The model 
shows the transition from the metallic state V < V c to a MH 
insulator for V > V c with V c = 2t. Our study is focused on 
systems with non-zero t', for which there is no exact solution. 
We choose such a system due to the lack of p-h symmetry, 
since for t' = the LSR is automatically satisfied for n = 1/2 
D8Q. In the following we study t' jt = 0.4, 0.2 in order to have 
substantial deviation from the p-h symmetry, but at the same 
time to maintain the simple momentum distribution for non- 
interacting fermions. 

Let us first consider the opening of the MH gap Ao at 
n = 1/2. It is well known [14] that within the t-V model the 
gap opens at V = 2t being exponentially small for V > 2t+ 
and nearly linear in V for V > At, see Fig. 1. We analyze the 
effect of t' on Ao by performing the exact diagonalization of 
chains with N = 14, 18, 22, 26, and 30 sites using the Lanc- 
zos technique. Chains (with periodic boundary conditions) 
were chosen to have odd number of electrons N e = N/2, 
since this leads to a non-degenerate ground state. The MH 
gap is then determined via 
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To obtain the gap in thermodynamic limit as Ao = 
limjv— >oo A (TV) we perform the finite-size scaling for N — 
14 - 30 using A(N) = a + b/N + c/N 2 . In the metallic 
regime, V < 2t, the gap scales expectedly with a ~ 0, and 
at large V ^S> 2t with 6^0. In general parameters a, b, c 
were obtained from least-square fit. Scaled gap Ao = a for 
t' = 0.4< is shown in Fig. [T] together with the exact Bethe 
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ansatz fl4ll result for t' = 0. The result reveals that the gap 
for non-symmetric case deviates only slightly from the exact 
result for symmetric case. 
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Figure 1: (Color online) Scaled gap Ao for non-symmetric case t' = 
OAt (points), and exact Bethe ansatz results for t' = (line). 



The (retarded) Green's function G(k, u>) at T = is defined 



as 



G(k,w) 



oo 

-I / dte^ + ^(0\{cic k (t)} + \0), 



where jj, is the chemical potential. We are studying finite sys- 
tems at fixed number of electrons, N e = N/2. Clearly, the 
position of \i within the MH gap is crucial for further discus- 
sion of LSR Si- Within this approach the correct choice is 
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For LSR ReG(fc,0) is important. In order to start with 
a general perspective we present in Fig. [2] the evolution 
of G(k,0) with increasing V. For non-interacting system 
ReG(fc,0) has a singularity at kp — 7r/2 persisting appar- 
ently within the metallic state, V < 2t. This is the normal 
Fermi-liquid (in ID Luttinger-liquid) behavior connected with 
the poles (singularities) of the spectral function A(k,u>) = 
—ImG(k,oj)/ir approaching ui = for k — > kp. With in- 
creasing V/t and the opening of the MH gap the behavior 
changes qualitatively. Within the gap A(k,u>) = 0, hence 
G(k, 0) has no singularity and goes through zero smoothly at 
k ~ 7r/2. Moreover, for V > At, G(k,0) becomes small for 
all k. This happens because at V 3> t, A(k,u>) consists of 
two nearly equal weights at approximately ±V and their con- 
tributions to real part of GF almost cancel each other at the 
chemical potential, therefore G(k, 0) ~ l/V 2 111 311 . 

Let us now focus on the LSR and its breakdown. The con- 
tent of LSR is the precise locus of k — kp, where G(fc,0) 
changes sign |U101 5 called Luttinger wave vectors 0|. Ac- 
cording to the LSR, for spinless model at n = 1/2 one should 
generally have kp — ir/2, if the topology of the electronic 
band is not changed qualitatively (which could happen, e.g. 
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Figure 2: (Color online) Evolution of ReG(fc, 0) with increasing in- 
teraction V for t' /t — 0.4 and for half-filled system of N — 26 sites 
with allowed wave vectors (points) and indications of kF ~ tt/2 
(dashed line) and (full line). Points are connected with straight 
lines. 



for t' > 0.5t). From Fig. [2] we note that kp is indeed near 
7r/2, however, even without finite-size scaling a small devia- 
tion kp it/2 may be observed for V > At. More accurate 
analysis with the finite-size scaling is presented below. 

To determine kp in the thermodynamic limit N — > oo, we 
perform the finite-size scaling of results at various N. Our 
procedure is the following. Due to periodic boundary con- 
ditions, allowed are k = 2-kI /N . Since we work with non- 
degenerate case with odd N e = N/2, k = tt/2 does not 
appear directly for any system. The closest are, however, 
k~(N) = tt/2 - n/N and k+(N) = tt/2 + n/N. Next we 
determine [i using Eq.© for given N. Then, the GF is evalu- 
ated via ED using the Lanczos algorithm for T = dynamical 
quantities IU7I1 , i.e., from Eq.fO, we have G = G c + G a , 

G c {k,uj) = (0|4(w + n + E° N - H)- 1 ^}, (4) 



and analogous for G a . 

GFatfc ± (7V), 



In particular, we calculate values of 



G ± (N) = G N (k ± (N),0), 



(5) 



where Gat(/c, ui) stands for the GF of a system-size N. These 
values are used to calculate kp in the following manner. For 
each N, we evaluate the mean value G(N) — (G + (N) + 
G~(N))/2 and the difference AG(N) = G+(N) - G~(N), 
which lateron serve for the evaluation of the derivative 
dG(k,0)/dk. Next, we perform scaling of both, G(N) and 
AG(N), to obtain their values in the limit N — > oo. 
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Finite-size scaling of G(N) is performed by assuming 



(6) 



Quadratic term Ci-4y is included in analogy with A(Af). For 
V ~ 4i dominant is the linear term 1 /AT, but with increas- 
ing V, parameter b\ in Eq. © decreases and ci term be- 
comes more important. In Fig. Rvalues of G(N) and obtained 
scaling from least square fit is shown for t'/t = 0.4 and for 
two values V/t — 4, 10. As it is seen from Fig. [3] the rele- 
vant limiting value is G = limjv— ►«> G(N) = ai. Note that 
ai = G(tt/2, 0) should be zero according to the LSR. Obvi- 
ously our finding in Fig. [3] that a% ^ is the indication that 
the LSR is violated. 
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Figure 3: (Color online) Scaling of G(N) for two values of V/t and 
t' = OAt. 



For the estimate of ki, also the scaling of AG(N) is needed. 
G(k, 0) within an insulator is a continuous function of k, 
hence AG(N) goes to zero as N — ► oo, so that the proper 
scaling function is AG (AT) = b 2 /N + C2/N 2 . Least-squares 
fits for the same parameters as in Fig. [3] are shown in Fig. 
[4] We notice, that for large V > lOt the linear term 6 2 /N is 
dominant, while for V ~ 4i quadratic correction with c 2 also 
becomes relevant. 

For further analysis the derivative of GF at ir/2 is relevant 
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To calculate the Luttinger momentum k^ for which GF 
changes sign we use linear approximation of GF near tt/2 as- 
suming that ki does not deviate appreciably from tt/2 being 
indeed the case. The result for kp may in our approximation 
be written as 



G(n/2,0) 



k =-- 

L 2 dG(k,0)/dk\ k=n/2 
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(8) 



Final results for two values t'/t = 0.2, 0.4, and within 
the whole range of parameters V/t are shown in Fig. [5] with 
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Figure 4: (Color online) Scaling of AG(N) for two values of V/t 
and t' = OAt. 
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Figure 5: (Color online) Calculated Luttinger momentum Ul vs. V/t 
for two values of t'/t. 



corresponding error bars. Presented values of fei in regime 
V > 3.6t are calculated as described above via Eq. (H). In a 
window 2 < V/t < 3.6 results are not shown since finite- 
chain calculations become unreliable due to limited wave- 
vector resolution and deviations from the simple finite-size 
scaling behavior. On the other hand, within the metallic 
regime V < 2t, G(k, 0) has singularity near tt/2. Therefore it 
makes more sense to apply instead a similar scaling analysis 
to the inverse values of GF, (ReG(k, 0)) _1 , and to locate in 
this way k^. Results obtained in this way are shown in Fig. 
[5] for regime V < 2t. Our estimate of error bar in Fig. [5] is 
given as the larger value obtained either from standard devi- 
ations of parameters ai and 62 or from the difference of the 
scaled value without taking into account the smallest system. 

From Fig.|5]it is evident that at large V/t > 6, k^ substan- 
tially deviates from the LSR prediction kp = tt/2. The devi- 
ation Afci = kh — tt/2 from LSR saturates at large V ^> t, 
whereby its value scales with the asymmetry given by t'/t. 
At the same time, with the decreasing V/t — > 2 and with 
vanishing of the MH gap Aq, also Afcf, appears to vanish. 
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In fact, in the regime 2 < V/t < 6 the LSR deviation Akp 
seems qualitatively to follow the variation of Ao. This is in ac- 
cordance with our observation that inverse derivative (— I/&2) 
has a similar behavior to Ao as a function of V. For large 
V > 6t, both 62 as well as 02 behave as oc l/V 2 , hence Ak^ 
approaches a constant value which could be evaluated via the 
method of moment expansion [13]. On the other hand, the 
analysis of data within the metallic phase for V < 2t does 
not show (within our accuracy) any deviation from the LSR 
kp = 7r/2. Hence, our results are consistent with previous 
confirmations of the LSR in the metallic phase away from half 

miing HEi. 

In conclusion, our results clearly show that the LSR is vio- 
lated in the Mott-Hubbard insulator within the ID generalized 
t-V model where the p-h symmetry is broken via the introduc- 
tion of the n.n.n. hopping t' ^ 0. Although we concentrated 
only on two values of t'/t — 0.2, 0.4, the behavior is quite 
generic whereby the violation Akp seems to scale with the p- 
h asymmetry, at least for modest parameter t'/t. It should be 
stressed that substantially larger t'/t can perturb qualitatively 
the topology of the noninteracting band and nk which makes 
the interpretation more difficult. An important finding is that 
the violation Akp scales as the MH gap Ao on approaching 
the metallic transition at V — V c , while within the gapless 
(metallic) phase we do not find any evidence for the LSR vio- 
lation. 

Discussing possible generality of our finding, we first stress 
several advantages of the extended ID t-V model analyzed 
above. Being a spinless model it allows for an accurate 
enough study of largest systems using the ED (compared, e.g., 
to the Hubbard model or multiband models), in our case up to 
N = 30 sites. Large span of sizes then allows for a reliable 
finite-size scaling and extrapolation to N — » 00. 

Another important ingredient is the absence of a phase with 
a long range order in the phase diagram at n ~ 1/2. The 
latter is due to the ID character of the model and supported 
by the exact solution at t' = M14I1 which is not strongly 
perturbed by moderate t'. Thus we are dealing solely with 
the MH metal-insulator transition at V = V c and with its in- 
fluence on the LSR. Note that in analogous D > 1 models 
studied so far, e.g., the 2D Hubbard model on a square lattice 
or triangular lattice 1 13| as well as more general MH 
insulators |@], the MH insulator is mostly accompanied with 
an onset of a long-range magnetic ordering, while doped MH 
insulators can show ferromagnetic order etc. An appearance 
of the long-range order clearly spoils the translational invari- 
ance and thus the validity of the LSR Ja] or at least requires 
the reformulation of the latter. 

In view of above discussion, the studied model is just the 
prototype example of a model with the MH transition and one 
can expect that in analogy the LSR would be generally vi- 
olated within the MH insulators without the p-h symmetry, 
even more if the MH transition is followed with an ordered 
state breaking the translational or some other symmetry. 



Quite open question is, however, the origin of the break- 
down of the LSR, or more precisely which part of the orig- 
inal proof JH |H oh becomes invalid within the MH insula- 
tor. The basic argument QJJ] invokes the existence of the func- 
tional constructed via the perturbation expansion in the inter- 
action strength. The required adiabatic connectivity 10, S @] 
to noninteracting fermions can be clearly questioned at the 
metal-insulator transition which at least represents a nonan- 
alytic point at n = 1/2 in the extended t-V model, even at 
t! = [14]. Moreover, it is evident that zeroes G(k, 0) = 
in the MH gap at the same time require divergent self energy 
E(/c, 0) — ► 00 JllQ] which enhances doubts on the existence 
of an appropriate functional. The question is important in 
connection with the emerging conserving approximations for 
strongly correlated systems based on a construction of such 
a functional but also on dynamical-mean-field methods based 
on approximations for S(fc, to) |2^ 21]. In any case, a deeper 
understanding of the limitations of the validity of LSR is still 
missing and can provide new insight in the physics of strongly 
correlated electrons. 
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